Continuum Coupling and Spectroscopic Properties of Nuclei 



o 
o 

(N 
oo 



> 

in 
in 
o 

o 
\o 
o 



X 



N. Michel,^' ^'^ W. Nazarewicz,^^^''^ M. Ploszajczak,^ and J. Rotureau^'^'^ 

^ Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996 
® Physics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831 
Joint Institute for Heavy Ion Research, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, Tennessee 37831 
^ Institute of Theoretical Physics, University of Warsaw, ul. Hoza 69, 00-681 Warsaw, Poland 
^ Grand Accelerateur National d'lons Lourds (GANIL), 
CEA/DSM - CNRS/IN2P3, BP 55027, F-14076 Caen Cedex, France 
(Dated: February 9, 2008) 

The nucleus is a correlated open quantum many-body system. The presence of states that are 
unbound to particle emission may have significant impact on spectroscopic properties of nuclei. In 
the framework of the continuum shell model in the complex momentum-plane, we discuss salient 
effects of the continuum coupling on the spectroscopic factors, one-neutron overlap integrals, and 
energies of excited states in ^He and ^*0. We demonstrate the presence of non-perturbative threshold 
effects in many-body wave functions of weakly bound and unbound states. 
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One of the main frontiers of nuclear science is the 
structure of short-lived, radioactive nuclei with extreme 
neutron-to-proton asymmetry. Such nuclei inhabit out- 
skirts of the chart of the nuclides, close to the particle 
drip lines, where the nuclear binding comes to an end. In 
this context, the major challenge for nuclear theory is to 
develop theories and algorithms that would allow to un- 
derstand properties of those exotic physical systems pos- 
sessing new and different properties [EQ- To this end, 
a unification of structure and reaction aspects of weakly 
bound or unbound nuclei, based on the open quantum 
system (OQS) formalism, is called for. 

The nuclear shell model (SM) is the cornerstone of our 
understanding of nuclei. In its standard realization 
SM assumes that the many-nucleon system is perfectly 
isolated from an external environment of scattering states 
and decay channels. The validity of such a closed quan- 
tum system (CQS) framework is sometimes justified by 
relatively high one-particle (neutron or proton) separa- 
tion energies in nuclei close to the valley of beta stabil- 
ity. However, weakly bound or unbound nuclear states 
cannot be treated in a CQS formalism. A consistent de- 
scription of the interplay between scattering states, res- 
onances, and bound states in the many-body wave func- 
tion requires an OQS formulation (see Refs. 0, Q and 
references quoted therein) . Properties of unbound states 
lying above the particle (or cluster) threshold directly 
impact the continuum structure. Coupling to the parti- 
cle continuum is also important for weakly bound states, 
such as halos. A classic example of a threshold effect 
is the Thomas-Ehrman shift ^ which manifests itself in 
the striking asymmetry in the energy spectra between 
mirror nuclei. Another example is the so-called helium 
anomaly [7|, i-e., the presence of higher one- and two- 
neutron emission thresholds in *He than in ^He. 

In this work, we investigate the impact of the 
non-resonant continuum on spectroscopic properties of 
weakly bound and unbound nuclei. In particular, the 
one-neutron overlap integrals and energies of excited 



states are studied using the SM in the coniplex fc-plane, 
the so-called Gamow Shell Model (GSM) By 
explicit many-body calculations that fully account for a 
coupling to scattering space, we demonstrate the pres- 
ence of a non-perturbative rearrangement in the wave 
function with a significantly low ang ular momentum 
single-particle (s.p.) component [TlLfl^. 

GSM is the multi-configurational SM with a complex- 
k s.p. basis given by the Berggren ensemble 13] con- 
sisting of Gamow states (poles of the s.p. S'-matrix; 
sometimes called Siegert or resonant states) and the non- 
resonant continuum of scattering states. For a given par- 
tial wave,(£, j), the scattering states are distributed along 
I ■ 

the contour L_^_ in the complex momentum plane. The 
Berggren basis is generated by a finite-depth potential, 
and the many-body SM states are the linear combina- 
tion of Slater determinants built from the s.p. Berggren 
states. The attractive feature of the GSM is that it can 
treat simultaneously both continuum effects and many- 
body correlations between nucleons due to the configu- 
ration mixing. All details of the formalism can be found 
in Ref. i. 

Single-nucleon overlap integrals and the associated 
spectroscopic factors (SFs) are basic ingredients of the 
theory of direct reactions (single-nucleon transfer, nu- 
cleon knockout, elastic break-up) (lj,il3- Experimen- 
tally, SFs can be deduced from measured cross sec- 
tions; they are useful measures of the configuration mix- 
ing in the many-body wave function. The associated 
reaction-theoretical analysis often reveals model- and 
probe-dependence 0, raising concerns about 

the accuracy of experimental determination of SFs. In 
our study we discuss the uncertainty in determining SFs 
due to the two assumptions commonly used in the stan- 
dard SM studies, namely (i) that a nuclcon is transferred 
to/from a specific s.p. orbit (corresponding to an ob- 
served s.p. state), and (ii) that the transfer to/from the 
continuum of non-resonant scattering states can be dis- 
regarded. 
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The SF for a given reaction channel is defined as a real 
part of an overlap between the initial and final states |l5j| : 

s^^{Wt-^®u,,Y-\^i-f, (1) 

where \^''/-^) and l^-;^-") are the wave functions of nuclei 
A and ^—1, and uij is the wave function of the trans- 
ferred nucleon (normalized to unity). The radial part 
of uij represents the one-nucleon radial overlap integral. 
Using a decomposition of the channel in the com- 

plete Berggren basis, one obtains: 

S'=n(^{^i^\\a+{bmit-^f^ , (2) 

where aj(6) is a creation operator of a s.p. basis state 
that is either a discrete Gamow state or a scattering 
continuum state. The normalization constant M is de- 
termined from the condition that when l^*^"*) — 

ll^l"!-!^ ]"''*)• Since Eq. ^ involves summation over 
all discrete Gamow states and integration over all scatter- 
ing states along the contour , the final result is inde- 
pendent of the s.p. basis assumed. This is in contrast to 
standard SM calculations where the model-dependence 
of SFs enters through the specific choice of a s.p. state 
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FIG. 1: The real part of the overlap {®He(g.s.)|['^He(g.s.) (g) 

^3/2]" )^ in three different SM calculations as a function 
of the energy of the 0p3/2 resonant state. The solid line 
(GSM) shows the full GSM result. The dotted line (HO-SM) 
corresponds to the SM calculation in the oscillator basis of 
OP3/2,Opi/2. The dashed line (GSM-p) shows the GSM result 
in the pole approximation. The imaginary part of is shown 
in the inset. 

In the following, we shall consider the ground state 
(g.s.) SF of ^He corresponding to the channel pHe(g.s.)® 
^3/2]*^ J where the single-neutron channel consists 
of the Op3/2 resonant state and the non-resonant 
scattering states on the complex contour ll^''^ . The 



s.p. basis is generated by a Woods-Saxon (WS) poten- 
tial with the "^He" parameter set |^ which reproduces 
the experimental energies and widths of known s.p. res- 
onances 3/2j~ and 1/2^ in ^He. The GSM Hamiltonian 
is a sum of the WS potential, representing the inert ''He 
core, and the two-body interaction among the valence 
neutrons. The latter is approximated by a finite-range 
surface Gaussian interaction 9] with the range ^=lfm 
and the coupling constants depending on the total angu- 
lar momentum J of the neutron pair: V^^'^ = -403 MeV 
fm^, Vq^^ = -392 MeV fm^. These constants are fitted 
to reproduce the g.s. binding energies of ^He and ^He m 
GSM. 

In order to investigate the continuum coupling, the 
depth of the WS potential is varied so that the 0^3/2 s.p. 
state (the lowest P3/2 pole of the S- matrix), which is also 
the g.s. of ^He in our model space, changes its character 
from bound to unbound. The valence space for neutrons 
consists of the 0^3/2 resonant state, complex-momentum 
P3/2 scattering continuum, and non-resonant pi/ 2 scatter- 
ing states along the real-fc axis. We do not consider the 
Opi/2 broad resonance explicitly in the basis as it plays 
a negligible role in the g.s. wave function of ^He 0. For 
both _L_|f -contours, the maximal value for k is fc,nax =3.27 
fm~^. The contours have been discretized with up to 60 
points and the attained precision on energies and widths 
is better than 0.1 keV. 

The calculated SF pHe(g.s.) 0^3/2]°^ in ^He is shown 
in Fig. ^as a function of the energy of the 0^3/2 pole. 
The SF strongly depends on the position of the pole: for 
^0pv2 < it decreases with eop.^^^ while it increases for 
^Op3/2 ^ 0- the one-neutron (In) emission threshold in 
^He, eop3^2=0, the SF exhibits a cusp. At this point, the 
derivative of SF becomes discontinuous, and the coupling 
matrix elements between the resonant OP3/2 state and the 
non-resonant continuum reaches its maximum Our 
calculations are consistent with the early estimate of the 
behavior of the cross section near the threshold energy by 
Wigner ll.]. The quickly varying component of SF below 
the In threshold in ^He behaves as (— e^^)^"^/^. Above 
the threshold, SF is complex; the real part behaves as 
(ej^)^"*"^/^ while the imaginary part, associated with the 
decaying nature of '^He, behaves as (e^f)^"^/^. 

To assess the role of the continuum, both resonant and 
non-resonant, the GSM results are compared in Fig. 
with the traditional SM approach (HO-SM) and with the 
GSM in the pole approximation (GSM-p). In the SM- 
HO variant, the harmonic oscillator (HO) basis contain- 
ing 0p3/2 and Opi/2 states with fuj — 41 A" MeV was 
used. The s.p. energies were taken from the GSM and 
their imaginary parts were ignored. The resulting SF 
varies little in the studied energy range and no thresh- 
old effect is seen. It is interesting to note that in the 
limit of an appreciable binding, the 0^3/2 wave function 
is fairly well localized, the importance of the continuum 
coupling is diminished, and the SM-HO result approaches 
the GSM limit. In the GSM-p variant, only the OP3/2 and 
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Opi/2 resonant states are retained in the basis, i.e., the 
non-resonant scattering states are ignored. Here, the g.s. 

of ^He is described by an almost pure [Op3/2 (8) 0^3/2]°^ 
configuration and the resulting SF is close to one in the 
whole energy region considered. A dramatic difference 
between GSM and GSM-p results illustrates the impact 
of the non-resonant continuum. 
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FIG. 3; Similar to Fig. Q except for the overlap in the ex- 
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*0: (^«O(0+)irO(3/2+)®d3/2l"^>'. The 
first derivative of the SF in the neighborhood of the eodg^a^O 
threshold is shown in the inset (a) while the inset (b) displays 
the imaginary part of . 
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FIG. 2: Solid line: one neutron radial overlap integral for 

{^He(g.s.)|[^He(g.s.) (gips/a]"^}^ in GSM for three energies (- 
0.8, 0, 0.8 MeV) of the OP3/2 resonant state. Dashed line: 
radial wave function of the OP3/2 resonant state of the WS 
Hamiltonian adjusted to reproduce the In separation energy 
of ®He in GSM. The relative difference between corresponding 
(r^) (in percents) is also given. 

Figure 121 displays the associated In radial overlap in- 
tegral calculated in the GSM for three energies of the 
Op3/2 resonant state. With eopg^., increasing, the overlap 
becomes more diffused. For comparison, we also show 
the radial wave function of the OP3/2 resonant state of 
the WS potential with a depth adjusted to reproduce 
the In separation energy of ^He calculated in GSM. The 
agreement between s.p. wave functions and many-body 
overlap integrals is excellent. Asymptotically, both quan- 
tities fall down exponentially with a decay constant de- 
termined by the In separation energy of ^He. The effect 
of the non-resonant continuum is seen in a slightly better 
localization of GSM overlaps. (The squared radii of over- 
lap integrals are reduced by several percent as compared 
to those of s.p. states.) 

To investigate the dependence of SFs on the orbital an- 
gular momentum i of the S'-matrix pole, we show in Fig. 
13 the SF for the excited 0^ state of ^^O in the channel 
[^^0(3/2^) (8) ^3/2]°^- Here, the s.p. basis is generated 
by a WS potential with the "I'i'O" parameter set that 
reproduces experimental energies and widths of the 5/2^ 
and 1/2'^ bound states and the 3/2^ resonance in ^^O. 
As a residual interaction, we took the surface delta in- 
teraction with the coupling constant Vo=-700 MeV 
fm'^, which was fitted to reproduce the g.s. binding en- 
ergy of ^^O relative to the ^^O core. The valence GSM 
space for neutrons consists of the 0^5/2, Is 1/2, and c?3/2 



Gamow states and the non-resonant complex ^3/2 con- 
tinuum. The maximal value for k on the contour L"^'^ 
is fcmax =1.5 fm"^. The contour has been discretized 
with 45 points and the resulting precision on energies 
and widths is ~ 0.1 keV. 

The behavior of SF shown in Fig. Ois similar to that 
of Fig. ^ except the variations are much smaller and the 
threshold behavior is different. Namely, the SF is contin- 
uous; it is its derivative that exhibits a cusp around the 
0^3/2 threshold. Again, this is consistent with the gen- 
eral expectation that for £=2 Re{S^) below the In thresh- 
old of ^^O behaves as (— Cj^)"^/^ while above the thresh- 
old i?e(52) {Im{S^)) should behave as (eji)^/^ {{ejef^^)- 
The associated In 1^3/2 radial overlap integrals (not dis- 
played) are extremely close to the s.p. resonant 0(1^/2 
wave functions; the relative difference in square radii is 
around 1%. 

Coupling to the non-resonant scattering continuum 
may change significantly energies of the many-body 
states close to the particle-emission threshold 0. We 
now shall examine the interplay between the continuum 
coupling and spin-orbit splitting in ^He and ^He using 
the same GSM Hamiltonian as before. Since we are in- 
terested in the first excited l/2~ state of ^He, an explicit 
inclusion of the Opi/2 Gamow state is necessary. Gon- 
sequently, the Berggren basis consists of the OP3/2 and 
Opi/2 resonant states and their respective complex scat- 
tering contours L^^^ and L'^^^. We take fcmax =3.27 
fm~^ for both contours, which were discretized with 27 
points. The GSM eigenvalues for ^He have been found 
using the Density Matrix Rcnormalization Group tech- 
nic|ue |T9l | adopted for the non-hermitian GSM problem 

The energy difference AE — E-^^^- ~ ^3/2^ '-'^ 
lowest states in ^He and ^He, calculated as a function of 
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FIG. 4: The energy difference, AE, of 1/2^ and 3/2^ states 
in ^He and ^He as a function of the spin-orbit strength. The 
gray line shows AE in ^He. In this case AE is equivalent 
to a s.p. spin-orbit splitting. The full GSM result for ^He 
is given by the black solid line, and dotted and dashed lines 
mark HO-SM and GSM-p results, respectively. Experimental 
splitting in ^He is indicated. 



the spin-orbit term in the one-body potential, is shown 
in Fig. ^ In ^He, these two states in our model space 
are the s.p. Gamow states Op3/2 and Opi/2, so that their 
energy difference is a spin-orbit splitting. In the absence 
of the residual interaction, the same splitting is obviously 
obtained for ^He. 

To investigate the importance of the non-resonant con- 
tinuum, we have performed the complete GSM calcula- 
tion in the full model space, as well as the GSM-p and 
HO-SM calculations. In the pole approximation, AE 
follows closely the behavior for ^He, albeit it is slightly 
larger. The calculated GSM-p states in ''He are almost 
the same as in the absence of the residual interaction, so 
the ^He and ""He curves differ only by a small shift due to 



the residual interaction. In the HO-SM variant, the resid- 
ual interaction gives rise to a configuration mixing, and 
AE in '^He is slightly reduced as compared to the original 
spin-orbit splitting in ^He. In the full GSM calculations, 
the coupling to the non-resonant continuum comes to the 
fore. The energy splitting in ''He varies steeply with Kso, 
contrary to other curves. The difference between GSM 
and GSM-p results reflects different correlation energies 
in the 1/2^ and 3/2 j~ states due to a coupling to the 
non-resonant continuum. This means that effects of the 
one-body spin-orbit interaction are effectively modified 
by the structure of loosely bound or unbound states, so 
that the same interaction can lead to a different value of 
AE in neighboring nuclei. 

In summary, our OQS calculations demonstrate the 
importance of the non-resonant continuum for the spec- 
troscopy of weakly bound nuclei. Firstly, the GSM study 
shows that the behavior of SFs around the particle emis- 
sion threshold exhibits a characteristic non-perturbative 
behavior; the predicted near-threshold reduction has, in 
our model, nothing to do with short-range correlations. 
Secondly, the continuum coupling strongly impacts the 
energetics of excited states (Thomas-Ehrman effect), an 
effect that goes well beyond the traditional CQS descrip- 
tion. Interestingly, the effect of the non-resonant contin- 
uum on In overlap integrals is minor: the single-particle 
approximation often used in SM studies seems to work 
very well. 
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